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^ . Abstract 

CN ■ Extensions of the MSSM could significantly alter its phenomenology at the LHC. We 

study the case in which the MSSM is extended by an additional U(l) gauge symmetry, 
which is spontaneously broken at a few TeV. The production cross-section of sleptons 
is enhanced over that of the MSSM by the process pp — > Z' — > ££*, so the discovery 
potential for sleptons is greatly increased. The flavor and charge information in the 
resulting decay, £ — ► £ + LSP, provides a useful handle on the identity of the LSP. With 
the help of the additional kinematical constraint of an on-shell Z' , we implement a 
novel method to measure all of the superpartner masses involved in this channel. For 
certain final states with two invisible particles, one can construct kinematic observables 
bounded above by parent particle masses. We demonstrate how output from one such 
observable, mnj can become input to a second, increasing the number of measurements 
one can make with a single decay chain. The method presented here represents a new 
^ class of observables which could have a much wider range of applicability. 



1 Introduction 

Supersymmetric extensions of the Standard Model, with M so f t ~ 1 TeV, are probably the 
most theoretically elegant solutions to stabilize the hierarchy between Mp and A wea k. The 
minimal supersymmetric standard model (MSSM) maintains the same gauge symmetries 
as the Standard Model, while introducing superpartners for the Standard Model particle 
content. In the past two decades, "standard" experimental signatures of low energy super- 
symmetry have been carefully studiedpQ. 

The MSSM is the minimal implementation of low energy supersymmetry. However, 
minimality should not be taken as a fundamental guiding principle in the search for new 
physics. In fact, the particle content of the SM advocates, if anything, non-minimal physics. 
In particular, we consider it reasonable to anticipate extensions of the gauge sector as well 
as the Higgs sector of the MSSM. If the gauge structure or matter content of the MSSM is 
extended, the phenomenology could be significantly different. There will be novel features 
deserving attention. New techniques and observables will have to be developed to extract 
the full information about the underlying model. 
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Before outlining the main results of our study, we briefly remark on the current status of 
the study of experimental signatures of the MSSM. As mentioned above, many now classic 
signatures have been studied pQ. Recently, special attention has been paid to a set of bench- 
mark models[2]. In particular, in SPSla, it has been shown [HI [4J that mass measurements 
at the LHC could be achieved to a very high accuracy 

While such studies could be instructive, we remark that the well-studied benchmarks 
should not be regarded as generic points in the MSSM parameter space. Consequently, 
the methods of measurements employed probably only have limited applicability and the 
conclusions could be misleadingly optimistic. In a recent study, [5], it has been shown that 
there are degeneracies associated with both discrete ambiguities (such as LSP identity) and 
larger uncertainties (such as slepton masses) in such measurements for a generic point of the 
MSSM parameter space. 

If supersymmetry is discovered at the LHC, one of the biggest challenges will be the study 
of the properties of the electroweak-inos and the sleptons. In the MSSM, the production 
of electroweak-inos is usually dominated by the cascade decays of color-charged particles. 
Such events will typically have a large number of jets, which makes the properties of the 
electroweak-ino difficult to study. As shown in [5], copious production of leptons in SUSY 
signals, typically associated with on-shell slepton production and decay, will greatly enhance 
our ability to study the properties of these superpartners and eliminating degeneracies. The 
direct production of sleptons pp — > Z*/j* — > ££* that decay to electroweak-inos suffers 
from a lower rate as well as a large Standard Model background [6j [7j. Although sleptons 
and electroweak-inos are easy to study in benchmark scenarios such as SPSla, where nig > 
Mi 2 > mi > M L sp, this is not true generically. 

In our study, we consider the LHC phenomenology of extensions of the MSSM, with 
special focus on its electroweak-ino and slepton sector. For concreteness, we focus on one 
typical possibility of such an extension, an extra U(l)', m z > ~ 0{1) TeV, which couples to 
both quark and lepton supermultiplets. Such an extension is fairly generic as it is present 
in many GUT/string motivated top-down constructions [8]- [18]. For most of our study, we 
will consider, as an example, [/(1)b-l- Being the unique non-anomalous global symmetry 
of the Standard Model with generation-independent charges, it is perhaps the most likely 
extension to the gauge sector. We will consider more general possibilities in the discussion 
of discovery reach. We will demonstrate that the channel pp — > Z' — > U* greatly enhances 
the discovery reach of £, and that copiously produced sleptons give an interesting handle on 
the identity of the LSP. Roughly, this only requires raz> > 1m~ t The result of this study is 
presented in section HJ 

Measuring the masses of the superpartners is usually quite difficult, as most of the kine- 
matical observables only measure their mass differences. The unknown momenta of neutral 
LSP's lead to undetermined kinematical variables, hindering reconstruction of the event. 
Guesses of unknown variables are usually unreliable since there are several of them. Such 
a difficulty is expected to persist in the pp — > Z*/"y* — > ££* channel, which has 3 unknown 
variables. The existence of an on-shell Z ! provides one additional kinematical constraint and 
should enable us to do better. One of the main results of this paper is the development 
of a new method which fully takes advantage of such a constraint. This allows us to com- 
pletely determine the slepton mass and the LSP mass with properly chosen observables. Our 
method is presented in section HI 
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For most of this study we will assume for simplicity that the sleptons are degenerate 
in mass. We remark that we should expect a certain amount of splitting between the left- 
handed and right-handed sleptons, at least from the effect of RGE running from a high 
scale. Such effects are more prominent if the overall mass scale of the sleptons is low. In 
models such as gauge mediation where there is a significant contribution from a Z' B _ L , [IH] 
we expect a larger universal contribution to both left-handed and right-handed sleptons. 
Nevertheless, we emphasize that the effect of left-right splitting is important and deserves 
further study. A detailed consideration of this issue is outside the scope of this paper. The 
methods we introduce for mass measurement could be modified, possibly by introducing 
additional observables, in order to deal with this further complication. In the conclusion 
of this paper, we will further argue that this effect does not impact the discovery reach as 
long as Z' decays into both left and right sleptons are still allowed. If one of them becomes 
heavier than m^/2 however, the signal significance should drop accordingly. 

In our study, events are generated at the matrix element level using COMPHEP-4.4.3 [2D] 
and piped through PYTHIA 6.3 [21] for initial state radiation and hadronization. PGS[2"2"] 
is used as detector simulation^ 

In section [2], we investigate the reach at the LHC for sleptons in models with an extended 
gauge sector and compare it to the MSSM for certain benchmark scenarios. We show how 
to use the Z' to determine the LSP identity in section [3l In section HI we discuss mass 
measurements for these benchmark scenarios. 



2 Discovery 

In this section and for most of this the paper we concentrate on the production channel 

pp -> z' -> 1 1 -> e + r + $ T (1) 

where I can be a sneutrino as well as a charged slepton. We also include the MSSM process 

pp -> Z*/-f* 1 1 -> £+ r + $ T (2) 

for comparison. 

In general, we consider an extra U(1)b- x l, which couples to the SM fermions through 

£ D g i> q B -xL l^Z^ ip, (3) 

where qs-xL is the charge of the SM fermions under U{1)b- x l- 

In the MSSM, the Z* /j* mediated slepton pair production cross section falls sharply with 
s > 4m|, and therefore with increasing slepton mass. On the other hand, for < mz>/2, 
the production through Z' resonance is almost independent of mj, up to a very mild phase- 
space factor, leading to a great enhancement in the discovery reach. In figure [TJ we display 
the improvement of the charged slepton (one flavor) production cross section over the MSSM 
for our benchmark scenario, with a Z' B _ L at 2.0 TeV and Qb-l — 0.25. 

1 The version of PGS used for this study also includes modifications made by S. Mrenna and J. Thaler 
for the LHC Olympics. 
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-- MSSM 




Figure 1: The cross section for one flavor of charged slepton pair production in the MSSM 
(dashed curve) and our benchmark scenario with a Z' B _ L at 2 TeV and gs-L — 0.25 (solid 
curve). 

In our analysis we preselect events with an opposite sign same flavor dilepton pair and 
use a jet veto (where jets are identified using a cone algorithm with R=0.7 and a pr cut of 
10 GeV). This significantly reduces SM background at high p T . In our background analysis, 
we include W/Z pair production generated by PYTHIA, and in the case that the final state 
lepton is an electron we consider the effect of W + j — > fake e + e~ (which we find to be 
negligible). The background analysis assumes \fN statistics, which allows us to find a rough 
estimate for the slepton reach at the 5cr level. We use data samples of 100 fb~ . 

In the preselected events we consider two observables, the $ T of the event and the 
of the softer lepton. As expected, we find that with these choices signal can generically be 
distinguished from background with relative ease by cutting on either of those observables. 
This is demonstrated in figure [2j We also include the MSSM case for comparison and see 
that it is difficult to separate from the background. Direct production in the MSSM was 
studied in detail in [6l [7] . 

2.1 LHC reach for sleptons in the MSSM 

In the MSSM, the available decay modes of sleptons depend on the identity of the LSP, 
and are especially sensitive to whether there is a chargino nearly degenerate in mass with 
the LSP. Therefore, we study the discovery reach in several different scenarios with different 
LSP's. We discuss briefly the reach in the MSSM for each scenario in this section. 

Our first scenario (I) has a bino LSP, with MSSM parameters given by Table [TJ Sleptons 
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Figure 2: Preselecting final state signatures of a same flavor opposite sign dilepton pair 
using a jet veto, comparison of the differential cross section for $ T (above) and the softer 
lepton px (below) in background, MSSM direct production and Z' resonant production in 
our benchmark scenario with m; = 200 GeV and vn^sp = 100 GeV. 
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are produced solely through pp — > Z*/Y — > £ £* and they decay to the LSP. Sneutrinos do 
not play a role in this scenario. For a single flavor (taken to be e) we find that the slepton 
discovery reach at 5a is < 300 GeV. 

Our next scenario (II) has a wino LSP. In this case, the left handed sleptons can decay 
to Xi thereby reducing the number of dilepton events compared to (I), while right handed 
sleptons can only decay through the effects of gaugino and higgsino mixing. The sneutrinos 
also produce dileptons as they decay to xf- We find that the 5a discovery reach is mj < 175 
GeV. The pattern of the leptonic signatures in this scenario is quite interesting, which we 
will discuss in detail in the next section. 

We also include a scenario (III) with a higgsino LSP which will be of interest when we 
attempt to determine the identity of the LSP in the following sections. 



name 


M x 


M 2 




tan (3 


I 


100 GeV 


2 TeV 


2 TeV 


10 


II 


2 TeV 


100 GeV 


2 TeV 


10 


III 


2 TeV 


2 TeV 


100 GeV 


10 



Table 1: The MSSM input parameters for the scenarios we consider. We decouple additional 
states by raising their mass above the squark and slepton masses. 

In order to demonstrate our signatures in cases with straightforward physics interpre- 
tations, we have chosen non-LSP electroweak-ino mass parameters so that the identity of 
the lightest electroweak-inos are pure gauge eigenstates and mixings induced by electroweak 
symmetry breaking are negligible. We have in fact studied scenarios where the electroweak- 
inos are split much more moderately, m non _LSP <^ 800 GeV, and found that our results in 
the next few sections are virtually unaffected. 

2.2 MSSM+[/(l) B _ a;i Scenarios 

Taking into account the mass and coupling bounds from LEP [231 121] and from the Tevatron 
[25] in U(1)b-xL scenarios, we will only be interested in models where ^y- > x (6 TeV). 

In a general study of U(1)b- x l, we consider the spectrum of scenario I, with sleptons at 
400 GeV and the Z' at 2.0 TeV, and vary both x and the coupling strength g. In figure [3] we 
display the 5a discovery contour at 100 fb _1 , our benchmark point and the exclusion contour 
from LEP data for this choice of Z' mass. At very small g, not enough Z's are produced to 
overcome background and at very small x the branching ratio of Z' —>■ ££ is too small so g 
has to be rather large for discovery. 

Next, we study the reach of slepton discovery as a function of the Z' mass, which is one 
of the most important factors determining the rate and hence the reach. For concreteness, 
we couple to scenarios I and II a U(1)b~l with g = 0.25 (and, naturally, x = 1). We scan 
over the mass of the Z' B _ L and find a great improvement in the reach for sleptons at 5a over 
the MSSM. Our results are displayed in figure HI At relatively low my, we find that sleptons 
can be discovered in most of the kinematically allowed region < mz>/2. Heavier mz> are 
rarely produced, and so one needs enough phase space to win over background. Therefore, 
the sleptons have to be light enough to be produced far from the kinematic threshold. 
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Figure 3: In a scenario with = 400 GeV and mz> = 2.0 TeV where the U(l) couples 
to baryon number with strength g and to lepton number with strength xg, one flavor of 
sleptons can be discovered at the 5er level at 100 fb _1 in the region to the upper right of the 
solid curve in the g — x plane. LEP data excludes anything to the upper right of the dashed 
curve. Our benchmark point is indicated by the solid dot. If we increase the mass of the Z', 
both the reach curve and the LEP constraint will shift upward and to the right. 




Figure 4: In bino LSP scenario and wino LSP scenario, we display the 5a slepton discovery 
reach at the LHC at 100 fb _1 in the presence of a U(1)b-l with g = 0.25 and x = 1 as a 
function of the Z' mass. For comparison, in the MSSM, the reach is < 300 GeV in bino 
LSP scenario and ra- t < 175 GeV in wino LSP scenario. 

3 Identity of the LSP 

In this section, we investigate the possibility of using the lepton information of the leptons 
from Z' decay to study the identity of the LSP. 

Depending on details of the mass spectrum, it is possible to study the identities of the light 
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Vn w = 2.28 x 10- 4 , Vn s = -0.9998 


wino LSP 


1.4 


7 x o ^ = 0.999 F x? 3 = -2.28 x 10" 4 
V x ± jW = -0.99996 U x ±-w = -0.998 


higgsino 
LSP 


3.3 


V y0W = 0.032 V y0§ = -0.018 

A,l ) A,X""^ 

= 0.024 = -0.013 
^ = -0.057 = -0.008 



Table 2: The correlation between the electroweak-ino mixing and the ratio of dilepton events 
to single lepton events in the three scenarios we consider with mi = 400 GeV for 100 fb~ 
of data. The ratio has been obtained for pure signal but includes detector effects. V x o^ etc. 
denote the relevant entries in the electroweak-ino mixing matrix. 

electroweak-inos by studying soft jets or the possible existence of displaced vertices in cascade 
decays involving nearly degenerate electroweak-ino states [261 [27]. The electroweak-ino can 
have a decay chain with large branching fraction into isolated, energetic leptons. These carry 
charge and flavor information and also provide a very powerful clean experimental probe to 
the identities of the electroweak-inos. In the MSSM, the best channel for such a study is 
X2 ~ > £ £ ~ > ^ + ^~ Xv However, as commented in the introduction, we should not regard 
this as a typical channel as it requires a certain arrangement of the spectrum. In our case, 
copiously produced sleptons from the Z' decay and £ — > t Xi will carry additional important 
information and provide a further handle on the LSP identity. 

In the following analysis we will be interested in the ratio of opposite sign dilepton events 
to single lepton events with a jet veto, restricting ourselves to one flavor (e ± ) as before. 
As we have already seen, in the case of a U(1)b-l, the dilepton signal dominates over 
SM background at high px- The single lepton signal has a larger background from singly 
produced W ± bosons but using a jet veto and concentrating on high px events still leaves 
us with a statistically significant excess, as we will show later in this section. 

The characteristics of the leptonic signature from the decay £ — > £ + LSP depend mostly 
on the bino and wino content of the LSP. A mostly bino LSP has the very distinctive feature 
that dilepton events greatly dominate over single lepton events, as charged sleptons always 
decay to \i y i a charged leptons. Any observed single lepton events with this process are due 
to detector effects. This can be seen in our example scenario I, as shown in Table [2j On the 
other hand, a wino LSP will offer the roughly comparable possibility of both dilepton and 
single lepton signatures. The slepton decaying directly into neutral wino LSP will produce 
only dilepton signature, just like the bino LSP case. On the other hand, there is a charged 
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wino state which is usually nearly degenerate with the LSP. Since it is very difficult to detect 
the existence of the process x — > X° + s °ft particles, we can effectively treat the chargino as 
the end of the visible decay chain. Slepton decaying processes £ + £~ — > t + x T + v + x° and 
z>z>* — > ^ + x T + v + X° wn l gi ye r i se to single lepton signatures, while vv* — ► t +X T 
can give rise to additional opposite sign dilepton signatures. Therefore, the ratio 

^ # of OS dilepton events 

^ # of 1 lepton events ' 

should give us a very clear handle distinguishing the wino and bino LSP shown in 

Table El 

The higgsino LSP case is more intricate. Since the Yukawa coupling of electrons is 
negligible, observables depend solely on bino/wino components of the lightest neutralino, as 
well as the wino component of the charged higgsinos, making the higgsino LSP scenario more 
difficult to distinguish from the other cases. In fact, we find the dilepton to single lepton 
ratio in our scenario III to be closer to the wino case (~ 1 : 1) than the bino case (> 100 : 1), 
as can be seen in table [2J If the gaugino/higgsino masses are changed, the amount of mixing 
will be affected and the signatures of scenario III can vary between the extreme cases of 
scenarios I and II. 

Having established an important difference between scenarios with different LSP identi- 
ties, we discuss the possibility of distinguishing these scenarios in the presence of background. 
For signal with single lepton, we use a background sample of diboson and single W ± produc- 
tion and look for events with a single very high-pr electron in the absence of jets. While we 
consider MSSM background as well, we find that after cuts its contribution is negligible. We 
study the single lepton production rate and find an excess over background of 4.2a in our 
scenario II and 2. 4a in scenario III for 100 fb^ 1 of integrated luminosity where the statistical 
significance of the excess lies in the region px > 500 GeV. Not surprisingly, scenario I does 
not give rise to any excess in single-lepton events so the higgsino LSP case lies roughly in 
the middle of the pure wino and bino cases. 

While the signal-only ratios in table [2] become much more uncertain due to background, 
we illustrate in figure [5] that scenarios I and II can still be distinguished due to the combined 
observability of single lepton and dilepton events, where using a jet veto we count single 
lepton events with p? > 500 GeV and dilepton events where the softer lepton has pt > 
400 GeV in 100 fb" 1 of data. With 100 fb" 1 of data, scenario III is not clearly distinguishable 
from scenario II. In fact, as we remarked earlier, the higgsino LSP case is expected to 
interpolate between the bino and wino LSP the gaugino content of the LSP varies 

between being bino-like and wino-like, while the coupling of the higgsino to electrons is 
irrelevant. 

While we have not taken advantage of them in our analysis, there are other potential 
differences between signal and background, such as the charge asymmetry in single lepton 
events from W ± production due to the pp-initial state which is absent in Z' initiated events. 
Exploring such effects fully can enhance the signal to background ratio in a more complete 
analysis. 

We mention here one additional example of the sensitivity of leptonic signatures to the 
structure of the electroweak-ino sector. One can have an e-/i non- universality coming purely 
from gaugino-higgsino mixing. In our scenario II, we find that for Jar, which does not couple 
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Figure 5: The number of dilepton [j)T2 > 400 GeV) and single lepton (pt > 500 GeV) events 
with a jet veto in 100 fb _1 of data. The background is given by the point close to the x-axis 
while the points away from the x-axis represent (from left to right) scenarios I, III and II 
respectively. The sizes of the error bars denote la using ^/N statistics, which include both 
the signal and the background. 

to the wino, the Yukawa coupling is just large enough for the branching ratio to Xi to be 
greatly enhanced through higgsino mixing while cr decays exclusively to Xi through its bino 
component. This is an interesting source of e/fi asymmetry that merits further study, but 
for us this just means that e + e~ is a better final state to consider than [i + [i~ . 

In summary, we see that the leptonic signature provides a very strong handle on the 
identity of the LSP. The higgsino LSP case or longer decay chains with more electroweak- 
inos could still give degeneracies but we expect a dramatic decrease in such degeneracies 
compared to the general MSSM. This agrees with the observation made in [5] in the case of 
on-shell sleptons. Since Z's generically give us a new source of on-shell sleptons, they help 
improve our ability to untangle the electroweak-ino sector. 

4 Measurements of Masses 

As argued above, in MSSM scenarios it is not generic for sleptons to be produced on shell 
copiously, and even if discovery is possible, one does not necessarily have enough statistics 
for mass measurements. In this section we will look in more detail into certain measurements 
made possible by the presence of a spontaneously broken U(l). 

4.1 General considerations 

In this section, we consider the case where the end of the SUSY decay chain is a stable 
electroweak-ino. For concreteness, we focus on the bino-LSP scenario. Due to the generic 
nature of the method we present here, we expect that it is straightforwardly applicable to 
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Figure 6: Momentum labeling conventions for Section HI 

other identities of the electroweak-ino LSPs. We will begin by reviewing the property of a 
set of generic p^-like observables, using rriT2 as an example. 

The rriT2 variable, as developed in [28], [29], [30l [31] , offers a straightforward way to calculate 
slepton mass in certain two-body decays. In our case, sleptons are pair-produced by the Z' 
with I — > I + x processes on either side of the decay as in figure By properly combining 
the lepton and missing energy information with the value of m Xl , one can measure m^. 

To construct m^, the unknown momenta q%, g 2 of the two x's are assigned in every 
way that satisfies the missing energy constraint qiT + q2T = p T - For each assignment of 
momenta, one constructs the transverse mass squared for both halves of the decay, 

^t(Pit, qiT, m Xl ) =m 2 e + m 2 xi + 2(E l T E* - p 1T ■ q 1T ) 

where Et = a/Pt + m2 ■ Since the transverse mass satisfies my < m^, the branch with 
a greater mf, gives a tighter constraint on m~ v Taking the greater of the two m^'s and 
minimizing this quantity over all possible assignments of q\ and q2 gives m<r 2 for the event. 
This is by construction less than the true transverse mass for one branch of the decay, which 
is in turn less than the mass of the parent particle, in our case m^. Thus, the distribution of 
m,T2 for all events has an endpoint at m^. 

The same information could be extracted from p? distributions of the visible particles. 
Very roughly, we expect the lepton p? distribution to peak near — m Xl . In principle, we 
could make this statement quantitative by simulating the decay process for various input 
masses and fitting to the resulting px distributions. The rriT2 endpoint does not carry more 
statistical weight than such a fit, but it has the practical advantage that it gives a quantitative 
measurement with a simple interpretation that does not require fitting to simulated data. 

To compute m T2 and reliably interpret its endpoint as the slepton mass, we must already 
know the LSP mass m xi . When m Xl is unknown, a free input mass M guess takes its place in 
the rriT2 equation: 

min [max {m% (p 1T , q iT ; M guess ) , m 2 T (p 2T , <\2T] M guess )}] (5) 

qir+q2T=p T 

We use a 130 fb _1 sample of selectron and smuon pair production (opposite sign same flavor 
dileptons, $ T , no jets or photons harder than 20 GeV) with M\ = 100 GeV and M~ t = 400 
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GeV. Knowing m Xl to be 100 GeV, and plotting the distribution of m^, we find the endpoint 
with a linear fit and measure = 405 GeV (see figure [7j). 
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Figure 7: m T2 for t -> i + \\ with 130 fb" 1 . Endpoint at 405 GeV. (m xi = 100 GeV, 
m~ e = 400 GeV) 

In the scenario where we know m Xl from another measurement, m^ 2 allows a quick and 
accurate determination of [28, 30J, even if we ignore the on shell Z' at the top of the 
decay. If m Xl is unknown, then the mT2 endpoint gives one constraint for the two unknown 
masses m Xl and mg, but it is not clear how to interpret it. The endpoint does not fall at 
either the slepton mass or the mass difference, and it does not vary linearly with the input 
mass, M guess . As we will show below, the presence of the Z' allows one to use rriT2 without 
this extra piece of knowledge. That is, we can measure both m f - and m Xl at the same time. 
To this end, we present a technique that determines m Xl to within 15 GeV. 

To understand how mT2 varies as a function of the unknown LSP mass M gucss , we gen- 
erated a sample of unreasonable luminosity, 7 ab _1 (100,000 events), and calculated m^2 for 
a wide range of input masses. One might hope that m T2 would change linearly as M guess is 
varied, but this only holds in the limit of large input mass, asymptoting to a line of slope 
1. This matches the behavior found in [3T] . Figure [S] shows how the accuracy in determin- 
ing mi depends on the uncertainty in m Xl . The discrepancy between at the correct 
value of M g uess and is a systematic effect due to radiation and detector resolution and is 
independent of luminosity. 

Thus, in the generic mr2 scenario, without outside of knowledge of m~ t or m xi , one obtains 
a curve in m^-m Xl space. With the extra constraint offered by a Z', one might hope to reduce 
the uncertainty in m Xl , collapsing the curve to a much smaller region. In this case, we have 

8 9i,2 ~ 2^ - lm LSPl =m LS p 2 - ln» fl =mj 2 ~ ^m z , ~ 1 m T2cndpoint = 2 Unknowns. (6) 

In a situation without a Z', we would have three unknowns, which we could take to be q^f 
and m xi . Given values for these unknowns, we could reconstruct the event up to a fourfold 
algebraic ambiguity from solving a quadratic equation for g x for each half of the event. In 
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Figure 8: ^T2 ondpoint for vs. M guoss for m Xl = 100 GeV. 100,000 events. 

the presence of the Z', however, we have a constraint on the total four-momentum. With a 
properly sensitive quantity, one could hope to show that only in a small region of m^-m Xl 
space does one sensibly reconstruct the Z' at its predetermined mass. 

In general, one could employ two strategies to achieve this. First, we could make some 
kinematical guesses about the unknowns and try to reconstruct the kinematics. As mentioned 
above, for any event there is the following hierarchy: mT2 < rriT < m£, where tut is the 
actual transverse mass of one branch of the event. For values of m^2 near the endpoint, one 
has approximately the correct my, so one might hope that the reconstructed LSP momenta 
are also approximately correct. This will be true provided that only a small region of LSP 
momenta is physically allowed after constraining = ^T2 cndpoint - Attempts to measure the 
bino mass with this technique, taking the 10% of events with rriT2 closest to the mri endpoint, 
were accurate to within only 100 GeV for 130 fb _1 . As shown in figure [9j the LSP momenta 
we reconstruct with this method only correlate roughly with the Monte Carlo truth, limiting 
the accuracy of this technique. As a second approach, we could take all possible values of 
the unknown variables and try to construct some observable which is bounded by the true 
values for the unknown variables. We will present a method based on this latter strategy. It 
works considerably better than the first, measuring the bino to within 15 GeV for 130 fb _1 . 

4.2 Constructing an Endpoint at mz> 

The basic approach is to use to compute as a function of M guess , and then to impose 
another constraint by demanding that the initial Z' be on shell. A step-by-step outline of 
our approach is as follows: 

1. Measure the Z' mass in an unrelated channel, such as Z' — > e + e~ . 

2. Guess Mg Uess , the mass of the LSP. 

3. Use m T2 (M gueS s) to compute the slepton mass. 

4. Compute the Z' mass. This is done by reconstructing every event in every possible way, 
picking the minimum allowed Z' mass for each event, then maximizing this minimum 
over all events. 
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Figure 9: Comparison of the actual p x of the LSP to the p x as reconstructed using the top 
tenth of events closest to the m T2 endpoint. 

5. Compare the Z' mass computed in step 4 to the actual Z' mass measured in step 1. If 
the answers are inconsistent, throw out this guess for the LSP mass M guess . 

6. Repeat this process for a range of LSP masses M guess . 

The definition of m^, eq.([5]), implies a "max-min" approach that can be applied very 
generally to processes with invisible final state particles: Given some unknown, M, construct 
an observable that is bounded above by the true value of M. Minimize the observable over 
all unknowns in each event, and plot the result for all events. The resulting distribution 
has an endpoint at the true value of m|§ This process can be applied sequentially, using 
the result of the first application as an input to the second. As an example, we first apply 
m T2 to compute m^(M guess ), then apply the strategy again to compute m z >- This gives us a 
measurement of one of the masses m x , m^/ if the other mass is known from another channel. 

The computation of rriT2 was described above. For given values of the unknown momenta 
qi,q 2 that are physically sensible (i.e. satisfy tjit < mi), we then use the slepton mass 
constraint (setting = ^T2 endpoiIlt ) to reconstruct my. Again, there is a fourfold ambiguity 
in rriz' resulting from the solutions to the two quadratic equations that determine q{ 2 . For 
each event, define the observable 



where the inner minimum is taken over the fourfold ambiguity, and the outer minimum is 
taken over all values of LSP momenta q\ and q 2 that reproduce the correct missing energy 

2 Clearly this method only works if the tail of the distribution is populated with events. Both tot2 and 
the max-min variable defined in this section satisfy this criteria, but extensions of m,T2 to processes with 
more than two invisible particles have fewer events near the endpoint 30J . Generally, minimizing over a 
larger number of unknowns makes it less likely for the endpoint to be populated. 




(7) 
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and obey < m T2 cndpoint - This observable clearly satisfies m™, m < mz>- Taken over many 
events, m™, in has an endpoint at the actual Z' mass (or, more accurately, at + Yz>) ■ 
Detector resolution, finite width of the Z', and the coarseness of our momentum sampling 
grid will smear the result. However, we nonetheless get an m^-like endpoint at the upper 
end of the Z' width. In the end, we are not actually interested in measuring mz', because it 
can be measured directly in another channel such as Z' — > £ + £~ . Instead, we use mz> as an 
additional constraint to determine the LSP mass m xi . 

We performed our analysis on an integrated luminosity of 130 fb _1 for M\ = 100 and 
Mi = 250 GeV. The Z' has mass mz> = 2 TeV and width Yz> = 27 GeV. The m™, m endpoints 
are shown in figure [10] for the correct input masses for M guess . We plot the value of the m™, m 
endpoint as a function of input mass for the two cases in figure [TTJ The uncertainty from the 
endpoint-fitting algorithm and Monte Carlo simulation (determined by repeating the analysis 
15 times with different random number seeds) is 27 GeV. There is an additional source of 
uncertainty we did not estimate, which is the expected position of the m™, in endpoint. In 
both scenarios we examined, this fell near the endpoint of the Z' width, + Yz> = 2.027 
TeV, but there is no reason to expect this to be exact. A determination of this uncertainty 
may loosen the bound we set, but we do not expect it to do so significantly. For both cases, 
over a range of a few hundred GeV, we sampled M guess at 10 GeV intervals. In the Mi = 100 
scenario, only for the correct guess mass of 100 GeV and 90 GeV did mz> + IV fall within 
the statistical uncertainty. For the Mi = 250 case, only the correct 250 GeV guess mass 
had mz> + IV in its uncertainty. These results are summarized in Table [3j In both cases 
we can determine the bino mass to within 15 GeV. An improved understanding of detector 
resolution effects could allow an even tighter bound. 

Finally, this can be compared to the original m<T2(M guess ) values to determine the slepton 
mass as well. For m Xl = 100 GeV we find = 405 with an uncertainty < ±10 GeV, and for 
m Xl = 250 GeV we find m~ t = 407 GeV with an uncertainty < ±15 GeV. These uncertainties 
include only the 27 GeV uncertainty mentioned above. We did not include the uncertainty 
in measuring the m<T2 endpoint. 

To determine the robustness of the analysis in the limit of low statistics, we compared 
uncertainties for M guess = 100 in the m xi = 100 case (figure [TTj) . Note that this uncertainty 
includes the systematic effect of the endpoint-fitting algorithm as well as the statistical 
uncertainty in the Monte Carlo. Down to ~ 250 fb~ , the analysis has the same sensitivity. 
It becomes progressively worse, but one can still constrain the Xi mass to a window 80-130 
GeV even for integrated luminosities below 40 fb _1 . Thus, one of the benefits of a max-min 
technique such as this is that such endpoints remain apparent even with only a few hundred 
events in the histogram. 



m Xl = 100 


M gucss = 80 


•Mguess 90 


M guess = 100 


m min Endpoint 


1998 


2015 


2028 


m Xl = 250 


M guess = 240 


M guess = 250 


M guess = 260 


m min Endpoint 


1988 


2026 


2055 



Table 3: Endpoints of the distributions of minimum m™ n calculated for each event (closest 
three to m z > + IV). All quantities in GeV. Uncertainty is ±27 GeV. 
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Figure 10: m™, m constructed in each event with 130 fb 1 (L: m Xl = 100 GeV. Endpoint at 
m^r = 2.028 TeV. R: m xi = 250 GeV. Endpoint at m™" = 2.026 TeV.) 




Figure 11: L: Endpoint of m™ n constructed in each event vs. M guess with integrated lumi- 
nosity of 130 fb -1 . R: Uncertainty in m™ n endpoint vs. integrated luminosity (m xi = 100 
GeV) 



4.3 Z' + Gravitino LSP 

So far, we have only considered the simple decay £ — > £ + x- Models dominated by cascade 
decays are more complicated, but the kinematics of the additional final state particles allow 
for a richer analysis. As an example, we consider a model with a massless gravitino and a 
heavy Z', where the slepton decays through a short lived bino NLSP, 

The events in figure [12] can be used to measure the masses of both the slepton and the NLSP. 
One approach is to treat the photons as "invisible" (by adding their transverse energy to ft?) 
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and proceed with the rriT2 analysis described above. However, this requires an artificially 
large number of events because it ignores the valuable kinematic information of the photons. 
Another approach is to devise a suitable max-min variable for this decay chain using the 
strategy described in the previous section. In order to demonstrate a different method, we 
instead use a weighting scheme based on the photon momenta to determine and m xi 
with a small number of events. 

1+ 7 




Figure 12: Decay of the Z' in a model with a massless gravitino and short-lived NLSP. 

Our analysis is similar to the measurement of the top quark mass in the dilepton channel 
at DO, tt — > jjuu£ + £^ [32J. In the case of the top quark, just 6 events were sufficient to 
measure the mass to an accuracy of 8%. In our case, there are two unknown masses instead 
of one, but there is also an additional constraint because each event starts with an on shell 
Z'. 

The spectrum, given in Table HI is motivated by gauge mediated SUSY breaking, but 
with the additional twist that the messenger sector is charged under {7(1)b_lI1 This gives 
a contribution to the scalar masses proportional to B — L, lifting above the heaviest 
electroweak gaugino. One of the telltale signs of gauge mediation is a kinematic edge in the 
dilepton invariant mass distribution from the decay x° ^£ + £^£ + x° + £, but this is 
forbidden in this model. On shell sleptons are produced only in the 

decay of the Z'. 



Particle 


Mass (GeV) 


Z' 
G 

Xi (mostly bino) 

£ 


2000 


100 
400 



Table 4: SUSY spectrum with a gravitino LSP used in this analysis. 

We studied an integrated luminosity of 100 fb -1 . Standard model background is negligible 
because of the two hard photons in the final state. Requiring two hard leptons, two hard 
photons, and no hard jets, we find 640 candidate events. 

3 This is similar to the "Harvard blackbox" model generated for the LHC Olympics [15] , 
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Given the mass of the Z 1 (which is easy to determine from another channel such as 
Z' — > e + e~), the kinematics of the event in figure [T2l are determined up to a single unknown. 
Therefore, for each such event, after imposing the constraints we are left with a one-parameter 
family of possible solutions for m f - and m Xl . Because the Z' has a finite width, and because 
the constraints are solved numerically on a coarse grid, the one parameter family of solutions 
is in practice a scattering of points in m f - — m xi space. 

All possible pairs (m ~, m l xi ) that solve the constraints are not equally likely. Ideally, each 
pair should be weighted by the probability that it would produce the observed event, 

W ideal ( {p obs } , mf , m xi ) = P ( {p obs }\m l = mf , m Xl = m Xl ) dp obs 

where {p bs} are the observed four momenta of the leptons and photons. The full probability 
function is difficult to calculate, so instead we use a simplified weighting function based only 
on the photon transverse momenta p\ and pf, 

W({p obs } } m}, m^) = P {p\,p 2 t I mi = m\, m xi = m l xi ) dp\dp\ . 

Photon momentum was chosen because the distribution of lepton p t from this decay is 
relatively flat. PYTHIA was used to generate the photon p t distributions for 45 reference 
models with 300 < m^~ < 600 and < m Xl < 200. Then, to compute the weighting function 
for each guess of (m~, my, the appropriate p t distribution was interpolated from the 45 
reference models. 

The total weight for each event is normalized to one. Finally, the weighted distributions 
for all 640 candidate events are added together, and the weighted frequency counts in each 
bin are interpreted as the "likelihood" of a given solution. The result is shown in figure [13j 
The maximal 5 GeV by 5 GeV bin is centered at m Xl = 107.5 GeV and m^r = 412.5 GeV, 
which should be compared to the Monte Carlo input masses m Xl = 100 GeV and m| = 400 
GeV. Clearly, this is a sensitive measurement of the slepton and NLSP masses; however, 
we have not done the Monte Carlo necessary to state reliable error bars. Even without a 
weighting function (i.e. W = 1 for all observables) , the maximum frequencies are found to 
be very close to the correct input masses. 

5 Conclusions 

Generically, we expect extensions of the gauge structure and the matter content of the 
MSSM. In this paper, we have studied the impact on supersymmetry phenomenology of U(l)' 
extensions of the minimal supersymmetric Standard Model. We have demonstrated that 
such an extension will give us a much stronger handle on the sleptons and electroweak-inos. 
Specifically, due to the enhanced slepton production cross section through pp Z' — > ££*, 
in comparison with the MSSM process pp — > Z* fj* — > ££*, we expect a greatly enhanced 
slepton discovery reach in this scenario. Moreover, with this additional source of on-shell 
sleptons, we have a much better handle on LSP identity. For example, simple signatures 
such as lepton counting reveal the existence of a chargino state degenerate with the LSP in 
the wino LSP scenario, distinguishing it from the bino-LSP case. 

With an additional resonance, a Z' in our case, we have more kinematical information 
to measure the masses of the superpartners involved in the decay chain. We developed a 
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Figure 13: Weighted frequency of m xi ,mi masses found by reconstructing 640 candidate 
events in every way consistent with kinematic constrains. The input masses are m xi = 100 
GeV and mj- = 400 GeV, and the maximal bin is centered at m xi = 107.5 GeV, = 412.5 
GeV. 

new method to take advantage of such a constraint. Using this type of max-min variable 
in the analysis of our model allows us to completely determine the masses of the sleptons 
and the LSP. This class of observables should have a much wider range of applicability in 
more complicated decay chains. Like standard edges and endpoints, m^Tike observables are 
particularly easy to implement because they do not require fitting parameters to a Monte 
Carlo simulation. 

In a generic model where new particles are produced in pairs and the final state has only 
two invisible particles, max-min variables can be devised that give additional constraints 
beyond the usual constraints provided by edges and endpoints. In the simplest case, pp — > 
££ — > Mxlxhi there are no standard edges or endpoints but the original m*r 2 variable gives 
one constraint. In a more complicated decay with more intermediate on shell particles, for 
instance the decay to gravitinos in Section I4~3"l the edge in the £7 invariant mass distribution 
gives one constraint on the two unknown masses and an appropriately designed max-min 
variable gives a second. As in the case of an on-shell Z', max-min variables can also be 
applied sequentially. The output of one variable can be used as the input to another applied 
higher in the decay chain. 

In our study we have made the assumption that the sleptons are degenerate for simplicity. 
Since the purpose of our paper is demonstrating the effect of a new set of signals and 
observables, rather than a comprehensive study, this assumption should simply be viewed as 
a useful first step. We do not expect small and moderate mass splittings of order 100 GeV to 
significantly affect the discovery reach since such a splitting should have little effect on the 
lepton pt and missing energy spectrum in figure [2J Making one of the sleptons heavier could 
in fact improve the discovery reach, as long as Z' — > £l,r are still both allowed, since the 
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slepton will decay into harder leptons. If one of the sleptons becomes too heavy for the Z' to 
decay into, we expect the signal significance to drop accordingly, depending on the change 
in the branching ratio to the remaining light slepton. For mass reconstruction, a method 
similar to the one we described, possibly by taking advantage of multiple end-points, should 
be applicable in the case of large L-R splitting. 

There are several other obvious directions to extend our study. First of all, one should 
also consider the case where the Z' decays into electroweak-inos. For example, this would 
be the case if the PQ symmetry were gauged and mixed with this U(l)'. Such channels will 
also allow us to have new windows into the structure of the electroweak-ino sector, which is 
generically difficult in the MSSM. 

We have studied the decay channel of the Z' into sleptons, but its decay into squarks 
is also interesting to consider. In this case, we expect to be able to extract additional 
information about the quark sector, complementary to that of the QCD production of such 
states. For example, if the Z' couplings to the quark states are chiral, we could have an 
additional handle on the left-right splitting of the squarks. 

Alternative extensions of the gauge structure of the MSSM will undoubtedly give rise to 
other novel features of phenomenology. It would be interesting to explore typical examples 
of such scenarios. 
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